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Abstract 

Many antimicrobial peptides (AMPs) have multiple antimicrobial immunity effects. One such class of peptides is temporins. 
Temporins are the smallest (AMPs) found in nature and are highly active against gram-positive bacteria. Nowadays, there 
was a rapid increase in the availability of the 3D structure of proteins in PDB (protein data bank). The conserved residues 
and 3D structural conformations of temporins (AMPs) were still unknown. The present study explores the sequence analysis, 
alpha-helical structural conformations of temporins. The sequence of temporins was deracinated from APD3 database, the 
three-dimensional structure was constructed by homology modeling studies. The sequence analysis results show that the 
conserved residues among the peptide sequences, the maximum of the sequences are 70% alike to each other. The second- 
ary structure prediction results revealed that 99% of temporin (AMPs) exhibited in alpha-helical form. The 3D structure 
speculated using RAMPAGE exposes the alpha-helical conformation in all temporins (AMPs). The phylogenetic analysis 
reveals the evolutionary relationships of temporins (AMPs), which are branched into seven clusters. As a result, we identi- 
fied a list of potential temporin AMPs which docked to the antiviral protein (MERS-CoV), it shows good protein-peptide 
binding. This computational approach may serve as a good model for the rationale design of temporin based antibiotics. 


Keywords Temporins - Alpha helical - Amino acid composition - Structure prediction - Phylogenetic tree - Protein-peptide 
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Introduction 
Antimicrobial Peptides 


Almost all multicellular organisms have AMPs to protect 
themselves from pathogenic microbes and acts as the first line 
of defence (Hancock and Chapple 1999). AMPs are an integral 
part of the body’s natural immune system which exhibits not 
only antimicrobial, anticancer, antitumor activity and but also 
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displays immunomodulating properties (Zasloff 2002). Major- 
ity of the mode of action studies performed using AMPs have 
cell membranes as their targets, and a variety of models have 
been proposed so far to explain the formation of membrane 
pores (Yeaman and Yount 2003). The disruptive membrane 
models such as “barrel-stave,” (Melo et al. 2009) “toroidal 
pore” (Hof et al. 2001) or “carpet models” (Steiner et al. 1988) 
provide a deeper insight into the interaction of amphipathic 
alpha-helical peptides with lipid bilayers, which are often 
referred as molecular electroporation (Matsuzaki 1999). The 
variations in cell membrane lipid composition of eukaryotic 
and prokaryotic organisms form as the primary targets for 
AMPs. The specificity of AMPs to the target cells relies upon 
the peptides favourable interaction with the microbial targets, 
which allow them to lyse a specific cell without affecting the 
host cells (Hancock et al. 2006). Net charge, amphipathicity, 
and hydrophobicity of AMPs also play a crucial role in the 
association of these peptides to target cellular membranes from 
exerting antimicrobial activity. Diverse antimicrobial peptides 
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have been reported from various ecological components such 
as plants, mammals, insects, and invertebrates. Currently, more 
than 2838 AMPs have been deposit in the antimicrobial pep- 
tide database. 


Temporins 


The temporins are small (10 to 14 amino acid residues), 
linear (AMPs) with a net cationic charge and an amidated 
C-terminus. These peptide classes are collectively known as 
temporins and have shown the highest activity among AMPs 
studied till date, against both human erythrocytes and bacte- 
rial strains. The mode of action of AMPs against bacterial 
cells has been extensively identified and documented; their 
antiviral functions were not fully understood (Skalickova 
et al. 2015). Recent findings have suggested the efficacy of 
various classes of AMPs and their specific mode of action 
against viral pathogens and the potential use of these AMPs 
as an antiviral therapy (Marcocci et al. 2018; Ahmed et al. 
2019). 

Due to the membranolytic effects of temporins, the patho- 
gens were difficult to develop resistance. Therefore, tem- 
porins are prospective aspirants for the design of imminent 
antiviral drugs with an innovative approach However, there 
is no structural information is available for the temporins 
family of AMPs. Therefore, we investigated the in silico 
structural information and protein-peptide interactions of 
temporins AMPs against MERS-CoV. These in silico find- 
ings are thought to help designing new types of peptide 
based antiviral drugs. 


Middle East Respiratory 
Syndrome-Coronavirus (MERS-CoV) 


In 2012, MERS-CoV was recognized in Saudi Arabia and 
belongs to the Coronaviridae family, mostly reported among 
the individuals of the Middle East (Milne-Price et al. 2014). 
MERS-CoV protein is a main target for therapeutic develop- 
ment to obstruct membrane fusion and viral entry (Du et al. 
2017). Since, there is no effective therapeutics for Middle 
East Respiratory Syndrome and this research attempts to 
find using peptide therapeutics by establishing a basis for 
the protein-peptide interactions through molecular docking 
analysis for the antiviral protein MERS-CoV. 


Materials and Methods 
Dataset 
Temporin family of (AMPs) was chosen for the present 


study. The amino acid sequence of (AMPs)belonging to 
the entire temporin family was retrieved from the APD3 
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antimicrobial database (https://aps.unmc.edu/AP/) (Wang 
et al. 2009, 2016; Wang and Wang 2004). A total of 121 
antimicrobial peptides of temporin family from various 
organisms were collected from the database. The sequence 
alignment was performed with the aid of MUSCLE pro- 
gram.(Edgar 2004; Zimmermann et al. 2018). Followed 
by performing Multiple sequence alignment (MSA), the 
temporin family was clustered based on sequence similar- 
ity. MMseqs2 and CD-Hit were used (https://weizhongli 
-lab.org/cdhit_suite/cgi-bin/index.cgi) for clustering of 
peptides (Steinegger and Séding 2017; Huang et al. 2010a, 
b; Li and Godzik 2006). Among the clustered peptides 
the conserved residues were determined using Weblogo 
(Crooks et al. 2004). The secondary structure prediction 
for the entire temporin family of peptides was performed 
to achieve the homology modeling. Further, the biofilm 
active temporin peptides were also predicted using the 
dPABBs module (Sharma et al. 2016). The half-life of the 
peptides was identified using the HLP, a server for predict- 
ing half-life of peptides (Sharma et al. 2014). 


Homology Modelling of Temporin AMPs 


The structural similarity of the sequences was performed 
using BLASTp search and the templates with similar 
structures in the PDB database were retrieved. Among 
the chosen 121 (AMPs) of temporins, 120 peptides were 
3D-modelled using Modeller software Version: 9.14 and 
PEP-FOLD for obtaining the three-dimensional structure 
of the peptides. The structure of one of the temporin pep- 
tide (PDB ID: 2MAA) 1s retrieved from protein databank 
structure (Sali and Blundell 1993; Fiser et al. 2000; Marti- 
Renom et al. 2000; Thevenet et al. 2012; Shen et al. 2014; 
Webb and Sali 2016). The modeled peptide structures 
were compared and evaluated with the standard peptide 
structure and validated using the Ramachandran plot at 
RAMPAGE server (Lovell et al. 2003). 


Phylogenetic Tree Construction 


A phylogenetic analysis based on the peptide sequence was 
carried out using the maximum-likelihood method with 
MEGA 6.0 (Tamura et al. 2013). The tree was evaluated 
using 1TOL (interactive Tree of Life) package (Letunic 
and Bork 2016). The pairwise distance was calculated and 
interactively visualizes their data in the form of the heat 
map (Babicki et al. 2016). 
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Molecular Docking 


To determine the molecular interaction of tem_a and high 
half-life (@¢n seconds) of temporin peptides with the antivi- 
ral protein (MERS-CoV) using Cluspro server, online based 
protein-protein, protein-peptide docking algorithm (Koza- 
kov et al. 2013, 2017; Vajda et al. 2017). MERS-CoV spike 
protein retrieved from protein data bank database (PDB ID: 
5X59). The selective temporin peptides were uploaded to 
the Cluspro server for docking analysis and the hydrogen 
bond interactions between the MERS-CoV antiviral protein 
with temporins are visualized by Discovery studio visual- 
izer software 2017(Dassault Systemes, BIOVIA Corp., San 
Diego, CA, USA). 


A Below 50% 


I) Between 51% - 60% 
Between 61% - 70% 
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Results and Discussions 
Sequence and its Physicochemical Properties 


The amino acid sequence of temporin AMPs are listed from 
various Rana’s family, and these AMPs are showing high 
hydrophobicity and positive net charge (Fig. 1). Temporin 
AMPs sequence fasta file supplementary |. The source and 
physicochemical properties of temporin AMPs were listed 
in Supplementary file 2. Most of the temporin AMPs are 
having more than 50% hydrophobicity (Fig. la), some of 
the temporin peptides (tem_hbl ‘FLPLLAGLAAKWFP’, 
tem_shl “FFFLSRIF’ and tem_lgd ‘FILPLIASFLSKFL’ ) 
has highest hydrophobicity percentage. In most cases, the 
hydrophobicity, length, net charges and other physico- 
chemical parameters determined the antimicrobial activity 
of a peptide. In specific the hydrophobicity and charge of 


B 


oOo-F- NW 


Fig. 1 Pie chart depicting distribution of hydrophobicity and net charge of temporin AMPs 
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Fig. 2 Conserved residues of temporin AMPs 
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<Fig.3 Multiple sequence alignment of temporin AMPs (Green colour 
indicates the sheet region; Magenta colour indicates the helix region; 
grey colour highlighted on temporins ID indicates the Bioactive pep- 
tides) (Color figure online) 


a peptide are critical factors determine the potency of an 
(AMPs) Any increase or decrease in the hydrophobicity 
index 1s directly proportional to the cationicity of a peptide. 
Most often the higher cationicity is attributing to greater 
binding efficiency, which results in effective binding of pep- 
tide over the bacterial surface (Malanovic and Lohner 2016). 
Once after binding the peptide takes its desired path to exert 
its action mechanism (Lee et al. 2002), while a decrease in 
hydrophobicity can reduce antimicrobial potency of the pep- 
tide (Kustanovich et al. 2002). During the synthesis of novel 
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Fig.4 The phylogenetic tree of temporin AMPs 


peptides, hydrophobicity is one of the perilous factors, for 
determining the range of target cells of AMP. Increasing the 
hydrophobicity of an AMP can alter the prospect of targets 
(Bahar and Ren 2013). 

The net charge of all temporin AMPs contains a posi- 
tive charge, and the average net charge is 2 (Fig. 1b). The 
selected AMPs of the net charge accounts for the ionizable 
groups of the peptide, resulting in the divergence from nega- 
tive to positive and vice versa, which is the foremost con- 
stituent for the initial interaction with negatively charged 
cell membranes. The average length of the temporin AMPs 
is 13 residues. The length of an AMP 1s also a factor that 
attributes to the bioactivity of a peptide. A minimum of 7-8 
amino acids are needed to form amphipathic structures with 
the hydrophobic and hydrophilic face on the opposite side of 
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Fig.5 A heat map showing the ism-67 
sequence identity between the 
temporin AMPs (red colour 
indicates the most similar 
pairwise distance green colour 
indicates the least similiar 
pairwise distance) (Color figure 
online) 
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a peptide molecule (Zelezetsky and Tossi 2006). The overall 
height of the amino acids in (Fig. 2) indicates the conserved 
residues at that position. We identified that non-polar amino 
acids like phenylalanine, proline, leucine, shows the con- 
served among the temporin AMPs. 

Sequence alignment of temporins AMPs represents the 
conserved residues among the peptide sequences, most of 
the temporin sequences are 70% similar to each other. The 
secondary structure prediction between the peptides shows 
the alpha helix and beta sheets, 98% of temporin AMPs 
exhibit alpha helical form. Some of the temporins display 
activity against the microbial biofilms, which are highlighted 
by grey colour shading (Fig. 3). From time to time, biofilms 
form a physical defense to protect bacteria from antibiot- 
ics. This might be due to slow development and reduced 
metabolic activity of a bacteria in such a population, the 
use of AMPs to inhibit biofilm formation could potentially 
be an efficient therapeutic option. In fact, due to the preva- 
lent mechanism of action of AMPs, that depends on their 
ability to permeabilize or form pores within plasma mem- 
branes, they also have a strong capacity to act on increas- 
ing bacteria.(Batoni et al. 2011).A recent study on temporin 
anti-microbial peptides demonstrated significant anti-HSV-1 
activity which was equivalent to that of Human cathelicidin 


G) Springer 


aie Fact 


rit -. a ell 


International Journal of Peptide Research and Therapeutics 


Values 


14.9806 
ff: 11.2355 


LL-37 and a temporin synthetic analogue [K3]SHa. In addi- 
tion, the study had emphasized its direct-acting mechanism 
for elimination of HSV-1(Roy et al. 2019). In an another 
in vitro anti-HSV-1 study using temporin B (TB), signifi- 
cant antiviral activity was observed which had resulted in a 
5-log reduction of the viral titer. The disruption of the viral 
envelope was also confirmed upon observation using a trans- 
mission electron microscopy. Further, the TB was found to 
disturb the HSV-1 life cycle in different phases, including 
the attachment, entry into the host cell, and post-infection 
stage (Marcocci et al. 2018). 


Homology Modeling of Temporin AMPs 
and its Validation 


Homology models of temporin AMPs have been constructed 
successfully for the identification of structural information. 
The result of all temporin peptides shows the a-helical con- 
formation in modelled structures, which is an important 
class of AMPs holding higher amphipathicity. Moreover, 
the a-helical conformation allows effective interaction with 
biological membranes (Huang et al. 2010a, b). The modelled 
3D structures of 120 temporins (AMPs) evaluated using 
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Fig.6 A scatter plot showing the half-life of the temporin AMPs (Green colour represent the minimum half-life, canal blue colour represent the 
neutral half-life, navy colour represent the maximum half-life) (Color figure online) 


RAMPAGE helps to interpret the allowing and disallowing Evolutionary Relationship and Recognition 

regions, which declares the evaluation of peptide structures of Half-Life of Tem porin AMPs 

in Supplementary file 3. 
We have constructed a phylogeny tree of all the 121 tem- 
porin AMPs (Fig. 4) to evaluate the evolutionary relation- 
ship between sequences. Maximum likelihood approaches 
are used to evaluate the Phylogenetic trees for a set of 
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Table 1 Energy of the docked complex based on cluster size (mem- 
ber), energy scores 


Peptide Target Representative Member Energy score 
Tem_a MERs-CoV _ Center 277 — 874.5 
Tem_Irnb 115 — 789.6 
Tem_Ip 223 — 925.0 
Tem_lIcsb 133 — 890.5 
Tem_cg3 133 — 861.0 
Tem_lec 191 — 1000.8 
Tem_If1 133 — 882.1 
Tem_eca 917 — 784.9 


species (Cho 2012), and it’s a robust statistical foundation 
that allows comparison of different trees, parameters and 
models. Based on the phylogenetic tree (Fig. 4), the tem- 
porin peptide sequences classified into five clades (CLADE 
A, CLADE B, CLADE C, CLADE D, and CLADE BE). 
Each clade node represents the individual colour to which 
the branched group has been recognized (Fig. 4). The Pair- 
wise distance matrix was calculated for all temporin peptides 
from MEGA software (Supplementary file 4) and a heatmap 
was constructed using HEATMAPPER, which shows the 
graphical representation of sequence pairwise between the 
temporin AMPs (Fig. 5). 

Plenty of the (AMPs) have entered clinical trials, but 
not many have been supported by the U.S. Food and Drug 
Administration to date due to the reason of short half-life. In 
this work, we predicted the short half-life of temporin AMPs 
by insilico method. (Fig. 6) shows the scatter visualization of 
predicted half-life values of temporin AMPs. Peptides like 
tem_Irnb, tem_lp, tem_lcsb, tem_cg3, tem_If1, tem_lec, 
tem_eca are showing the higher half-life values in seconds. 
The lists of the half-life of all temporin AMPs Supplemen- 
tary file 4. 


Evaluation of MERS CoV - Temporin AMPs 
and its Interaction Analysis 


The molecular docking of selective temporin AMPs with 
antiviral protein (MERS-CoV) was performed by using 
cluspro online server. For further evauation, we selected 
eight antimicrobial peptides, namely, tem_a, tem_Irnb, 
tem_lp, tem_lIcsb, tem_cg3, tem_If1, tem_lec and tem_ 
eca belonging to family of temporins. The modeled, vali- 
dated structure of tem_Irnb, tem_lp, tem_lIcsb, tem_cg3, 
tem_If1, tem_lec and tem_ecahave been chosen as ligand 
for this protein-peptide docking studies. The docked pose 
with the best fit was selected for each peptide-protein com- 
plex based on same binding region (Fig. 7). The cluspro 


result of MERS CoV — Temporin AMPs complex details 
presented in Table 1. 

Further, the protein-peptide complexes were visualized 
by Discovery studio visualizer, to display hydrogen bond 
interaction involved in the protein-peptide docking. The 
tem_a peptide shows good binding with the MERS-CoV, 
it forms six hydrogen bond interaction with GLU793, 
SER1095, LYS1021, ARG841 (Fig. 5a); tem_Irnb forms 
hydrogen bond interaction with GLN1084, ALA1080, 
THR1076 (Fig. 5b); tem_lcsb shows null hydrogen bond 
interaction with the target (Fig. 5c); tem_cg3 forms six 
hydrogen bond interaction with ARG841, GLN1097, 
LYS1021, ASN1028 (Fig. 5d); tem_lec shows 2 hydro- 
gen bond interaction with ARG841, LYS1021(Fig. 5E); 
tem_eca forms one hydrogen bond with SERI091 
(Fig. 5f); tem_lp forms three hydrogen bond interaction 
with SER1091, GLN1097 (Fig. 5g); tem_If1 shows good 
interactions, it forms five hydrogen bonds with LEU1094, 
SER845, ARG841, ASN1028 (Fig. 5h). Recently, Mustafa 
et al. 2019 had docked, investigated and validated the 
interactions between the various peptides (Cyclotides, 
Chicken beta defensin, Dermaseptin-s9, Human alpha 
defensin, Plectasin, Rat defensin, Rat alpha defensin) 
against MERS-CoV and has demonstrated prominant anti- 
viral activity. 


Conclusion 


In this present study, we have interpreted the sequence, 
modelled 3D structures, evolutionary relationship and 
docking studies with antiviral protein (MERS-CoV) of 
temporin AMPs through a homology modeling, phylog- 
eny tree construction and protein peptide docking. The 
3D structural of temporin antimicrobial peptides success- 
fully predicted and results shows the alpha-helical struc- 
tural diversity among the peptides. Phylogenetic analysis 
exposes that all temporin peptides are 70% evolutionary 
similar to each other. Some of the temporin peptides like 
tem_a and high scored half-life temporin peptides were 
docked to MERS-CoV protein, and results shows good 
interaction with the antiviral protein. These observations 
reinforce the idea that temporin AMPs may serve as a 
new therapeutic application for the future development. 
Furthermore, the findings reported above are preliminary 
and will certainly require experimental evaluation through 
in vitro and in vivo studies, which are crucial for their 
development into therapeutic targets against MERS-CoV. 
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